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1. Introduction 

Recently BaBar collaboration measured the semileptonic branching fractions B — )• Dtv and 
B ^ D*xv that are above their Standard Model predictions [1]. The experiment reports 

^^^^ = ^^{B^d17^) = 0.440 ±0.058stat. ±0.042,yst. , d-D 

for the decay with pseudoscalar D and the result was normalized with respect to the decay with light 
lepton in the final state in order to cancel Vd, and form factor parameterization in the theoretical 
prediction of this observable. Before this result was published a prediction based on the lattice- 
calculated form factors had been made [2] , 

/?(D)sM = 0.296 ±0.016. (1.2) 

The discrepancy immediately raised interest in the flavor community to explain it within one of the 
well-motivated NP models. On the one hand, the SM charged-current contribution to this decay 
hints that possible NP contributions should be present at tree-level. On the other hand, the SM 
prediction requires the S — )• D form factors whose knowledge from the lattice is limited to high-<7^ 
region where phase space is small. 

In this work we have revisited the theoretical prediction of B — ^ Dtv in the SM in a manner 
that maximally employs the available experimental information on the form factors from previ- 
ously measured B — Div where i stands for e,jLt [3]. Next, we parameterize beyond the Standard 
Model contributions to this decay mode in the effective Hamiltonian language, focusing on inter- 
actions that preserve lepton flavor universality and induce b ^ c transitions via scalar and tensor 
interactions. Assuming presence of either scalar either tensor operator we derive bounds on their 
respective Wilson coefficients at 1- and 2-a level, taking into account experimental and theoretical 
uncertainties. 



2. Differential decay width 

In the SM the amplitude for hadronic transition D — ^ P is given in terms of vector and scalar 
form factors, F+{q^) and Fo{q-), defined as 

{D{p')\cy,.b\B{p)) = (p, +p' - ""^^q,) F^{q'-) + ""^^q.Foiq'-) . (2.1) 
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The momentum transfer is denoted hy q = p — p' and when squared coincides with the invariant 
mass of the leptons, q^ = [pv + pf)^. In the SM, the differential width of B — ^ D£v decay valid for 
finite lepton mass m( is 



d^{B^D£Vi) _ Gl\Vcb\ 



dq-^ 



Tgo 



192713 



\Vcb\^^o\F+iq- 



2\\2 



cW)+ci{q^) 



(2.2) 



2 



B — ^ Dtv with minimal theory input 



Nejc Kosnik 



53^ 



+ 



1500 



1000 



500 




2 4 6 8 10 12 



12500 
10000 
3 7500 
5000 
2500 



0, 







4 6 



10 12 



Figure 1: Vector and scalar form factor weight functions for decay B DTVr are shown in the left plot. 
Dashed curve corresponds to CQ{q^), while the solid curve corresponds to c\{q^). In the right plot the same 
o(^^) plotted together with c^{q^) (red curve), the weight function in the case of massless muon in 
the final state. 



where 
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X{q'^,ml,mj)) = [q^ - {niB + muf][q^ - {niB - niDf] . 



The overall factor is defined by Eq. (2.2): 
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(2.3) 



(2.4) 



The scalar form factor Fq is multiplied by Cq(^7^) that is nonzero only when the lepton in the final 
state is massive. Therefore only the B — Dtv mode is sensitive to exchanges of charged Higgs in 
the context of two Higgs doublet models or alternative scenarios that induce scalar operators [4, 5, 
6, 7]. 

The ^^-dependence of the functions c^(^^) and CQ{q^) are shown in Fig. 1. On the left-hand 
side plot one observes that the contribution of scalar form factor is enhanced by large Cq(9^) that 
partially compensates smallness of FQ{q^). 



3. B ^ D£v with minimal theory input 

In this section we demonstrate, how one can maximally employ available experimental data 
to test for consistency of the SM in observable R{D). Semileptonic decays with light lepton {£ = 
e,li) have enabled extraction of \Vcb\- Usual approach in the literature has been to make the SM 
theoretical prediction of the ^^-spectrum at the maximal lepton recoil point q^^^ = [rtiB — mo)^, 
where the two hadrons are both at rest, and where the corrections to the heavy quark limit have 
been calculated. On the other hand, the experimental data close to q^^^ have virtually useless 
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Figure 2: In the left plot we show the binned spectrum of B — > Djiv below 8GeV^ (right vertical grey 
line) as measured by the BaBar collaboration [9]. Above 8GeV^ lattice QCD simulations in quenched 
approximatian (empty symbols) [10], and those in which A'f = 2 + 1 dynamical flavors are included (filled 
symbols) [11] are indicated. First grey line on the left denote the T threshold at = m^. Ratio of Z? — > D 
form factors obtained in the same lattice simulations are shown on the right plot. 



statistics due to small phase space, and experiments rely on the CLN shape of the form factor 
F^{q^) to fit the differential decay spectrum to data at low to moderate and extrapolate it to 

^max [8>9]. 

In order to maximally use the experimental data we propose a different path for both hght and 
heavy leptons in the final state of B — D£v. The vector form factor can be extracted from the BaBar 
measurement of the spectrum with light leptons in the final state. This procedure is viable below 
q^ = 8GeV^, where the errors in bins are small [9]. Prediction of B — )• Dzv in q^ G [SGeV^,^^^^] 
region will rely on the lattice QCD results for the vector form factor [10, 11] as seen in the left 
panel of Fig. 2. Notice that below the above procedure relies on the input value of Vet only at high 
q^ where lattice data is used. 

Now we can express both branching fraction entering R{D) as integrals over the three kine- 
matical regions denoted in Fig. 2 



!GeV^ 



^{B ^ DAtV^) = ^0 I c^4q')\V,bF+ 

Jml 

+|Kfe|'^o / "">^(g2)|f+.iatt(9')|W, (3.1) 

JSGeV- 



where q^^-^ = [ms — mj))^. The phase space integral for decay with heavy T contains the contribu- 
tion of the scalar form factor and is also spht at = 8 GeV^ 



!GeV- 



J ) I exp 



cW)+cl{q^) 



?max 
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dq 



dq^ . (3.2) 



We emphasize again that the partial decay rates at low-<7^ do not rely on the input value of Vcb, since 
\VcbF^{q^)\ is available experimentally. A suitable observable with this pleasing property should 
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be defined as a ratio of partial decay rates, i.e., 



ml 



^<8GeV 



q-<8 GeY- 



(3.3) 



and would depend only on |Fo/F+| tliat is well under control theoretically as we show below. 
The theoretical prediction of R{D), as defined in Eq. (1.1), requires in addition the SM value of 
Vcb and calculation of the vector form factor at high q^. For the latter we take the lattice results 
of Refs. [10, 11] and fit them to a dipole parameterization. The global fits of the SM to flavor 
observables yield \Vcb\ = 0.0411(16), a value that we use aboveSGeV^ [12]. 

The ratio of scalar-to-vector form factor is constrained to be 1 at = by construction and 
indicates, together with high results of lattice QCD, a linear behaviour (see Fig. 2) 

Fo{q^' 



F+iq' 



l-a q 



(3.4) 



We use the value a = 0.020(1 )GeV ^ that is consistent with different theoretical approaches (see 
e.g. [10, 11, 13, 14, 15]). The above procedure gives finally 



/?(D) =0.31 ±0.02, 
which is less than 2(7 below the BaBar result (1.1). 



(3.5) 



4. Constraint on NP effective Hamiltonian 



The SM effective Hamiltonian is extended to include operators that are scalars, tensors, or 
vectors, their dimensionless couplings labelled as gs, gr, and gy, respectively. We assume that the 
lepton flavor universality is respected by all couplings on the Lagrangian level. We do not invoke 
operators that contain a right-handed neutrino. 

J^eff = -VIGfV,!, [{cr^b){£LfVL)+gv{cr^b){£Lr^VL) 

+gs{li){cb){lRVL) +gT{pi){ca^vb){lRa^'''vL)\ +h.c. (4.1) 

All couplings scale as gv,sj °^ ni^/m^^, and /mnp is the new physics scale. The differential decay 
width in the presence of NP operators is 



dq^ 



\ycb\ 



5^0 \f+K 



,2m2 



+ 



CTviq^) Re 
l+gv)- 



l+gv|V+(^2) + |gr(M)l'4(9') 



FT{q\^) 



(l+^v)^*r(At) 
9^ gs{^) 



membilJ.) - mdlJ.) 



F+{q') 
'ci{q^ 



Foiq^ 



F+{q^ 



(4.2) 
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Figure 3: Regions of allowed values for gs{mb) and grimh), compatible with experimentally measured 
R{D). The small region within the solid, dashed and dot-dashed white curves correspond to the respective 
1-, 2- and 3-C7 compatibility with both ^('^''P) (B DrVt) and ^('"'P' {B-^DpiV^). The thick dot represents 
the Standard Model, namely gs.T{mh) — 0. 



The tensor form factor Fj{q^,iJL) is defined as 

{D{p')\ca^,b\B{p)) = -i {p^p'.-p'^py) (44) 

While the vector nature of weak currents has been thouroughly tested and are compatible with 
gv = 0, scalar g5(/x) and tensor g7-(/x) operators are not as constrained. ^5(/i) / induces the left- 
right operator and lifts the helicity suppression in B — Dtv and B — Dptv decays. Same operator 
is enhanced by the factor mj^/m^ with respect to the left- left (SM) contribution to the D° — 
mixing amplitude. Furthermore a noticable effect could also be seen in D — Vj decays that are 
governed by loops containing the down-type quarks and are therefore sensitive to 7^ [16]. 

Using R{D) alone we get a very loose constraint on gs{mh) while we require gy = gri^ib) = 
(contours on the left-hand plot in Fig. 3). Requiring in addition the compatibility of the theoretical 
expression for ^E{B — )• DjiVp^) obtained by using Eq. (4.2) and the measured value [9], restricts the 
allowed gs{mh) to a small region also indicated in Fig. 3. For example, when gs{mb) is real then 
the la compatibility with experiment allows —0.37 < gs{mi,) < —0.05, while the 3a compatibility 
amounts to —0.53 < gs{mb) < +0.20. The statistical error in Eq. (1.1) is treated as Gaussian, while 
the systematic errors and uncertainties with respect to the form factors are treated as uniform. 

If we allow for gr(wifo) 7^ in Eq. (4.2) then the possible values that are compatible with 
R{D) are those in the contour plot shown in the right-hand plot of Fig. 3. The needed tensor 
form factor has not been computed, to our knowledge, on the lattice nor in the QCD sum rules. 
The computation in the model of ref. [13] shows that FT{q^) /F^{q^) = 1.03(1) is a constant, in 
agreement with naive expectations based on the pole dominance. As before, R{D) alone is not 
constraining strongly the possible values of gj {mi,) whereas taking into account the constraint from 
measured ^{B — )• Dpiv^) [9] shrinks the allowed region, as shown in the right-hand plot in Fig. 3. 
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If Im grimb) = 0, we obtain 0.3 < grimh) < 1.5 and —0.6 < gTipih) < 2.1, from the requirement 
of respective 1- and 3a compatibility with both experimental R{D) and SS{B — )• DjXVp). 

5. Conclusions 

The result for R{D) could be an indication of new physics should the significance of incompat- 
ibility with the Standard Model raise to at least 3a. The compatibility with the Standard Model can 
be tested experimentally, with a minimal hadronic input, as discussed in this letter. Here we used 
the lattice QCD results for F^{q^) at larger because the full branching fractions were reported 
in Ref. [1]. Thus we have found R{D) = 0.31 ±0.02 where the significance of the discrepancy with 
Eq. (1.1) to be below 2a. 

If, instead of comparing the full branching fractions of both decay modes, the experimenters 
cut at about 8 GeV^, then the shape of the needed vector form factor including the factor of 
\Vcb\ could be reconstructed from the differential branching fraction of B — )• Dptv [4, 5]. The only 
theoretical hadronic quantity needed then is the slope of the form factor ratio (3.4), which is quite 
accurately known from lattice QCD with the values that agree with quaik models and with recent 
QCD sum rule studies. By using the vector form factor multiplied by \Vcb\ data from Ref. [9] only, 
and by integrating the decay rates up to 17^^ = 8 GeV^, we obtain 



[B^Dixv^] 



= 0.20 ±0.02. (5.1) 

r<8 GeV- 



AUowing for departures from the Standard Model, while keeping lepton flavor universality 
which has been experimentally verified to a very good accuracy [17], the measured R{D) and 
^{B — Dptv^) impose quite strong constraints on the new physics scalar and tensor effective 
couplings gsjinih). 
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